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Abstract 

This paper gives out the general solutions of variable coefficients Linear ODE and Riccati equation by way of integral 
series S{X) and TiX). Such kinds of integral series are the generalized form of exponential function, and keep the 
properties of convergent and reversible. 
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1. Introduction 

It is a classical problem to solve the n-th order Linear ODE : 



-u + ai{x)-j-^u + fl2(x)-^-^M + • ■ • + a„(x)u = /(x) 
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As we all known, 

1. if {fl„(x)} are all constants, Eq.(l) could be solved by method of eigenvalue ( Euler), or by exponential function 
in matrix form 

U = e^ ' ■ C + r^-' 



F (s) ds 



where C is a n X I constant matrix . 
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2. if {fl„(x)} are some variable coefficients, such as some special functions [Wang, P337,206] 

d^y 1 dy ^ ^ 

dx^ X dx x^ 

T d^y dy 

(1 - x2)-4 - 2x-^ + n(n + l)y = 

dx^ dx 

special function theory answers them. 



(Bessel Equation) 
(Legendre Equation) 



But when it comes to the general circumstances, the existing methods meet difficulties in dealing with Eq.(2) , because 
of the variable coefficients. In order to overcome it, two functions are invited : 



1.1. Definition 



r^x r^x r^t r^x r^t r^s 

£[X(x)]=/+ I X(t)dt+ I X(t) I X(s)dsdt+ I X(t) I X(s) I X(^)d^dsdti 
Jo Jo Jo Jo Jo Jo 

T{X(x)\^I+ \ X{t)dt+ \ I X(s)dsX(t)dt+ \ \ \ X{^)d^X(s)dsX{t)dt ■ 
V Jo Jo Jo Jo Jo Jo 



(4) 



It will be seen that such definition is reasonable and necessary. Clearly, when X{x) and J^X{t)dt are exchangeable, 
then 

S[X{x)] = e/o'^('>'" = T[X(x)] 

Besides, &{X) and 'FiX) extend some main properties of the exponential functions, such as convergent , reversible 
and determinant (see Theorem 3.1). In addition, a n x m matrix A{x) - {aijix)),^^ is bounded and integral in [0,b] 
means that all its element fl,j(x) are bounded and integral in [0,b]. 

2. Main Results 

Theorem 2.1. the general solution of the Linear ODE (2) is: 



U = S[A(x)] ■ C + &[Aix)] ■ [ T[-A(s)]-Fis)ds 

Jo 



(5) 



where C is a n X \ constant matrix . 

Theorem 2.2. For the bounded and integrable matrix , A{x) — {aij)„n, B{x) — (bij),„i„, P(x) — (pij)mn, Q{x) — {qij)n 
in [0,b], the general solution ofRiccati equation 



IS 



where 



or the other equivalent form: 
where 



— W + WPW + WB-AW - Q^Q 
dx 
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3. Solutions of Linear ODE 



3.1. Properties of&iX) and TiX) 
From the Definition(4), it holds that 



dx 

^ dx^ 



X{x) 
X{x) 



X-B 



X(x) 



r 



X{x) 
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Now, we will see more explicit properties of &(X) and T{X). 

Tlieorem 3.1 (Properties of&{X) and TiX)). IfX{x) is bounded and integrable, it holds that 

1. &{X) and 'Fix) are convergent; 
2. 

det£(X) = detr(X) = dete.t'^W'" = j:"-^^'^^, ^ ^,r£x(.,)d, 



3. S{X) and FiX) are reversible, and 



T(X)&{-X) = &{-X)r(X) = / 



(11) 



(12) 
(13) 



Proof. 1. Firstly, £(A) is convergent, since {fli(x)}"^j are bounded in [0,b]: 

3M > 0, s.t. lflt(x)| < M, Vx€[0,fe], k^\,2,---,n 



So 
(a) 



(b) 



(c) 



J A (/) t/f|| = max I flit (f) dt 
I Jo 



< M\x\ 



\^A{t) ^^A{s)dsdt'^ = maxW,i £ak{t) ^^ai{s)dsdt\ < !o^dsdt\ < jj{M\x\f 



A(t) \ A(s) I A(^)d^dsdt\\ — max \y I flit (0 I fli (*) I (^) c/f c/sfiff 
Jo Jo tt' Jo Jo Jo 



JO JO Jo 

r^t r^s 



l-J 



< n^M^ 



fif 

Jo Jo Jo 



d^dsdt 



< -^{Mixir 



(d) 

It follows that 



||£(A)|| < 1 + - 
n 



1 9 1 1 

nM\x\ + —{nMx) + —{nMx) + 



= 1 + -e' 
n 



,nM\x\ 



Clearly, S(A) is convergent. 
Similarly, FiX) is also convergent. 
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2. Vn X n matrix A{x), if trA{x) is bounded and integral , then 

which is a special case of Abel's formula[Chen] : IfWandB are n x n matrixes , s.t. 



then. 



—W^BW 

dx 



dstW^e" 



(14) 

(15) 
(16) 



Here we just take 2x2 matrix for verification: 









yi,i yi,2 




Let Y{x) = e 


A(x) 






,A(x)^ 








y2,i y2,2 





<^2,1 fl2,2 



SO 4- y = A ■ Fmeans that 

dx 



it follows 
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y2,2 





a\\ a\2 

fl2,l fl2,2 



yi.i yu2 
y2,i y2,2 



(17) 



— (dety) = det 

dx 



— det 



+ det 



y\.\ yi,2 

-J2,l -J2.2 



y2,\ y2,2 

ai.iyi.i + ai,23'2,i ai,i3'i,2 + ai,23'2,2 

y2,\ y2,2 



det 



)'l,2 



a2.\y\,\ + a2,2y2,l a2,l3'l,2 + fl2,2y2,2 



Al l det 



y\,\ 
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, y2,\ 


y2,2 



detY - trA ■ detY 



Thus, Abel's formula holds and &{A{x)) is reversible. 
By the times: 



detnX) = e/o'"-^W^' = (18) 
so, all we need to proof is 

deter^W'" (19) 

Because eii no longer satisfies Abel's formula (one reason is X and J^X(t)dt are unnecessarily exchange- 
able ) , we seek the other approach: 

V« X n matrix A, 3n x « reversible matrix P , s.t. 

P-^AP = diag{JuJ2, ■■■ , JA := J 
J is A 's Jordan matrix, 7,- is the Jordan block with eigenvalue Ai{x). 
It follows that 
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(20) 



So, 

p-ig'if = = = diag[e^\e^\ ■ ■ ■ ,e^'} 

Therefore 

dete'* = dete^ = = e"'-^ 

which yields 

dete/o'^*"'" = e/o'"'^(')^' 

3. Notice that V« x « matrix A, there exists a companion matrix A*,s.t. 

A-A''^A*-A^detA-I (21) 

so, i/det A , A /i invert ibie. 

Therefore, S{X) and ;F(X) are invertible. 
Furthermore, it holds that 

riX)6{-X) = £i-X)riX) = / (22) 

Because; 
(a) 



d_ 

dx 



T{X)8{-X) 



= ^r(X) ■ Si-X) + nX) ■ ^&{-X) = TiX)X ■ &{-X) - nX) ■ X3(-X) = 
ax ax 



So, 

nX)S(-X) = const. = [r(X)6{-X)]l^^^ = / 
(b) Due to the special property(21) of matrix, Eq.(22) is obtained. 

□ 

3.2. Proof of Theorem!. 1 

Proof. According to Definition(4) and Theorem 3.1 , it follows 

' ^&[A{x)] = A{x) ■ &[A{x)] 
ax 

d 

—G(x) = A(x) ■ G(x) + F 
dx 

where ^ 

G(x) = £[A(x)] • r r[-A{s)\-F{s)ds 
Jo 

-^G(x) =-^£[A(x)] ■ f r[-A(s)]-F{s)ds + S[A(x)]-r[-A{x)]-F{x) 
dx dx Jo 

=A(x) ■ £[A(x)] ■ r r[-A{s)]-F{s)ds + F 
Jo 

Clearly f/(x) = 6[A(x)] ■ C + &[A{x)] ■ £T[- A{s)] ■ F (s) ds is convergent. 

Moreover, since £(A) is reversible, U{x) is the general solution of Eq.(2). 
□ 



(23) 



because 
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Theorem 3.2. Assume that A{x) — {aij)„xn, B(x) — (bij)„,xm, P(x) — ipij)nxm ore bounded and integrable matrixes , 
and U{x) is the desired n Xm matrix. The Linear ODE : 



d 

— U^A(x)U+UB(x) + P(x) (24) 
ax 



has general solutions 



U{x) = £(A) 
where C is n x m constant matrix. 
Proof. Let U = £(A) ■ W ■ T(B), then 



Jo 



A(t))P{t)&{ - B{t))dt + C 



T{B) (25) 



So Eq.(24) could be reduced to 



— f/ = A(,x) U + UBix) + £(A) — W ■ TiB) (26) 
dx dx 

&(A)-^W -TiB)^ P (27) 
dx 

^ n-A) ■ P ■ Gi-B) (28) 

dx 

It's obviously that 

Wix)= f r[-Ait)]Pit)-S[-Bit)]dt + C (29) 
Jo 

C is n Xm constant matrix . □ 



or, 



4. Solutions of Riccati equation 

In mathematical investigation of the dynamics of a system, the introduction of a nonlinearity always leads to some 
form of the Riccati equation [ Watkins] : 



dx 



y + a(x)y' + b{x)y + c{x) - 



(30) 



But it is usually the case that not even one solution of the Riccati equation is known. In the following text, we try to 
give out solutions of Riccati equation in matrix form: 



dx 



W + WPW + WB - AW - Q ^ 



(31) 



where A(x) = {aij)„„, B(x) = {bij),nm, P{x) = {pij)mn, Q{x) = (^,7)™ . 

4.1. Proof of Theorem.2.2 

I'roof. 1. Firstly , define[Polyanin, Ch 0.1.4] 



W2 :^£(PW + B) 



(32) 



so W2 is reversible, if PW + B is bounded; 
meanwhile. 



Secondly, let Wi := WW2, so 



— W2^ (PW + B)W2 
dx 



(33) 



£Wi =£w-W2 + W-£W2 = £W-W2 + W- 



PW + B 



W2 = 



4-W+WPW+WB 



(34) 



so, with Eq.(31) and Definition (32), it holds 



— Wi = AWi +QW2 
ax 



(35) 



Take the relationship (33) and (35) into consideration. 



d 






A Q 




' Wi 


dx 






. P B 







(36) 



we can solve W\ and Wo- 

On the other hand, according to Definition (32) , it's obviously that 

W2|,=o = &{PW + B)|,=o = / 

so it goes without saying that 

mL=0 = (W2)|,=o = H'U=0 



(37) 
(38) 



We immediately obtain 



Therefore W - W\ ■ W^'^ '^^ the solution of Eq.(31). 
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^ W |,=o ' 
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_ W2 . 




P B 




/ 



2. Similarly, we can get 



-B P 
Q -A 



so, W = C/j ' ■ is also the solution of Eq.(31). 
3. But the two solutions are equivalence! That is. 



or 



t/2 • IVi - t/i ■ W2 = 



Because, according to Eq.(36) and Eq.(40) 



(39) 



(40) 



(41) 
(42) 



-W2 



UiP - U2A -Wi + U2 ■ AWi + QW2 - U2Q - UiB ■W2 - Ui ■ PWi + BW2 = 



(43) 



As a result, 

U2-Wi-Ui-W2 = const. = [U2 -Wi-Ui- W2]\^^^ = (44) 
which implied that two solutions are equivalence. 
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4. Uniqueness. If Eq.(31) has more than one solution,such as X(x),Y(x), under the same initial condition,i.e. 
X{Q) = Y{Q). Let W(x) = X{x) - Y{x). So it is clear that what we need to prove is equitant to show 



( d 

— W + WPW + WB-AW^Q 
\ dx (45) 

I W\,=o = 

has uniqueness solution W{x) = 0. 

Take advantage the proof steps we have established: according to step(39) and (35), 
Any solution ofEq.(45), such as W{x), it is reasonable to define 

W2 = S{PW + B), Wi = W ■ W2 

It follows that W2 is bounded , 

^Wi^AWi (46) 
dx 

and 

Wi = £[A] ■ Wi L,=o = £[A] ■ W|,=o = (47) 
Therefore, W{x) = Wi • W,"' = 



4.2. Simplify solutions ofRiccati equation by particular solution 

In the research of Riccati equation, particular solution plays crucial important role. Too much of works have been 
done. The first important result in the analysis of the Riccati equation is that if one solution is known then a whole 
family of solutions can be found [Watkins]. 



Theorem 4.1. The same conditions as theorem 2.2, Riccati equation 



— W + WPW + WB-AW-Q^Q (48) 
dx 



has the unique solution 



W + &{A-YP)- (WL,=o) 

where 



f 

Jo 



R{t){t)dt-(W |,=o) 



■T{ -[B + PY]) (49) 



Y is solution ofEq.(48) when W\x=o = 0, i.e. Y\x=o = 
R:=ri-[B + PY]) PBiA- YP) 



(50) 



Proof. 1. According to Theorem 2.2 , Eq.(48) has solutions. Take any one of it, such as Y, and let 

V^W-Y (51) 

It follows that 

VPV = (W - Y)P{W - F) = {WPW - YPY) - {W - Y)PY - YP{W - 7) = {WPW - YPY) - VPY - YPV 
£?.^8) (^^_d_y^ _WB + AW + Q]- [-^J -YB+AY+ Q]) - VPY - YPV 
= ( - j-^V + AV - VB) - VPY - YPV = --f^V + {A - YP)V - V{B + PY) 

(52) 

That is, 

d 

— V + VPV + V(B + PY)-(A-YP)V^O (53) 
dx 
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2. Obviously, &(a - YP) and T[-[B + PY]) are reversible , we may let 

V ^6[A-YP)-U ■r{-[B + PY]) 
Now Eq.(53) could be transformed into 



(54) 



S[a - YP) - —U ■T[-[B + PY]) + (A - YP)V - V(B + PY) 



+VPV + V(B + PY) - (A - YP)V = (55) 



or. 



3. Let 



— U + U 
dx 



r{-[B + PY]) ■ps[a- yp) 



U ^0 



R :^T[-[B + PY]) P-SiA- YP) 
According to Theorem 2.2, U has solution 



where 



= £ ■ =U+ dt] 

\ R Q I I \ Jo R / 



U |,=o 
/ 



U ^Wi- Wo' 




/ 



U |.v=0 

1 + £R{t)it)dt-U\,=o 



(56) 

(57) 
(58) 

(59) 



Now, Let's consider how to choose Y , so that both W and U | v=o are as simple as possible. It's clear that 

when y|,=o = 0, C/U=o = F|.v=o = W|,=o 

In this case. 



U = W|,=o • 



/+ R{t)it)dt-{WUi) 



(60) 



It should be noticed that 



/ + j^R(t){t)dt-(W\,=o) 



is reversible, otherwise 



Jo 



/+ R(t)it)dt-(W\,=o) = 



(61) 



which is clearly impossible. 



According to transformation(54), the solution of Eq.(48) is 



R{t) (t) dt ■ (W L-=o) ■T{ -[B + PY]) 





(62) 



where Y(x) = , if and only ifQ(x) = 0. 
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